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Abstract 

We present simulations within the framework of scalar-tensor the- 
ories, in the Newtonian limit, to investigate the influence of massive 
scalar fields on the dynamics of the collision of two equal spherical 
clouds. We employ a SPH code modified to include the scalar field to 
simulate two initially non-rotating protogalaxies that approach each 
other, and as a result of the tidal interaction, intrinsic angular mo- 
mentum is generated. We have obtained sufficient large values of J/M 
to suggest that intrinsic angular momentum can be the result of tidal 
interactions. 
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1 Introduction 



In recent years there has been achieved important progress in understanding 
the dynamics that led to the formation of galaxies. Two and three dimen- 
sional N-body simulations of galaxies and protogalaxies have been computed 
using up to a few millions of particles, giving a more realistic view of how 
galaxies, quasars, and black holes could have formed (Barnes & Hernquist 
1992, Barnes 1998). 

The Universe's composition at the time galaxies formed could be, theoret- 
ically, very varied, including baryonic visible and dark matter, non baryonic 
dark matter, neutrinos, and many cosmological relics stemming from sym- 
metry breaking processes predicted by high energy physics (Kolb & Turner 
1990). All these particles, if present, should have played a role in the structure 
formation. Then, galaxies are expected to possess dark matter components 
and, in accordance with the rotational curves of stars and gas around the 
centres of spirals, they are in the form of halos (Ostriker & Peebles 1973) 
and must contribute to at least 3 to 10 times the mass of the visible matter 
(Kolb & Turner 1990). 

Whatever the Universe composition was, protogalaxies were originated 
due to a spectrum of scale-invariant perturbations (Harrison 1970; Zel'dovich 
1972) that was present within the cosmological background at the beginning 
of structure formation; the inflationary cosmology is the most convincing 
scenario that explains its origin (Mukhanov et al 1992). Protogalactic struc- 
tures began to acquire some momenta, e.g. tidal torques (Fall & Efstathiou 
1980), because of local gravitational instabilities to provoke plenty of colli- 
sions, mergings, fly-bys, etc, between these original, cosmic structures. As 
a result of their evolution, galaxies, as we presently know them, must have 
formed. The dynamics of protogalaxies has been studied intensively, for a 
review see Barnes & Hernquist (1992). There has been much interest to un- 
derstand how galaxies acquired their present features, especially how their 
internal angular momentum (spin) has been gotten, J/M ~ O(10 30 ) cm 2 /s. 
A very important issue about it is how the transfer of angular momentum 
between protogalaxies took place to give rise to the observed elliptical and 
spiral galaxies with their known mass and rotational properties. As an ini- 
tial condition can be thought that protogalaxies were gravitationally isolated. 
However, there are some indications that the orbital angular momentum in 
spiral galaxies in pairs is few times larger than their spins, so pairs seem to be 
not dynamically isolated (Oosterloo 1993). Part of this angular momentum 
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could had its origin in the cosmic expansion (Caimmi (1989,1990), Andriani 
& Caimmi (1994)), where it has been computed the torques at the beginning 
of strong decoupling from the Hubble flow of spherical-symmetric density 
perturbations. Moreover, observations of spin angular momentum of various 
thousands of disc galaxies are compatible with the mechanism of generation 
of spin via tidal torques (Sugai & lye 1995). Other theoretical and numerical 
studies of evolution of angular momentum of protogalaxies from tidal torques 
are in line with observations (Chernin 1993, Catelan & Theuns 1996). 

In the present work we investigate how the transfer of orbital to spin 
angular momentum is achieved when two equal, spherical clouds pass by, and 
in some cases when they collide; this type of interactions are to be expected in 
the tidal torques scheme. Studies of interacting spherical systems very related 
to ours have been done using the Newtonian theory of gravity, including a 
number of topics: mergings (White 1978, 1979), mixing processes (White 
1980), simulation of sinking satellites (White 1983a), and mass and energy 
lost in tidal interactions (Aguilar & White 1985), among others. However, 
we made our calculations within the framework of scalar-tensor theories, in 
the Newtonian limit, to investigate the influence of masive scalar fields on 
the dynamics. 

This paper is organized as follows: In section 2 we present the Newtonian 
approximation of a typical scalar field theory. In section 3, we present our 
models of protogalaxies, the initial conditions, and the results. The conclu- 
sions are in section 4. 

2 Scalar Fields and the Newtonian Approxi- 
mation 

We consider a typical scalar field theory given by the following Lagrangian 

+ L M {g»v) (1) 
from which we get the gravity equations, 

R„-\^R = f (2) 



c = 



^ -0i? + ^(90) 2 -v(0) 

l07T 
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and the scalar field equation 



We expect to have nowadays small deviations of the scalar fields around 
the background defined here as ((f)) = 1. If we define = — 1 the Newtonian 
approximation gives (Helbig 1991) 

#oo = \v 2 h 00 = 4tt P - ^V 2 (4) 
V 2 0-m 2 = Snap (5) 

where we have done 



(j)V' - 2V 
3 + 2u 



m 2 <f) — m 2 k(j) 2 + . . . 



and a = 1/(3 + 2a;). 

The solution of these equations is 



= 2au\ 

h 00 = —2u — 2au\ (6) 



where 



u = 



E 



m a 



r-r 



ux = E^^exp[-|r-r a |/A] (7) 

a l r r a\ 

A = 1/m, where m is the mass given through the potential. This mass can 
have a variety of values depending on the particular particle physics model. 
The potential u is the Newtonian part and u\ is the dark matter contribution 
which is of Yukawa type. The total force on a particle of mass rrii is 

£ F = _I V /ioo = m,a • (8) 
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3 Protogalactic Cloud Models and Results 



Original protogalaxies could have very irregular forms, but we use as a first 
approximation spherical clouds for simplicity, and because this form seems to 
resemble the global shape of both visible and dark matter of many galaxies, 
i.e., taking into account their spherical halos (Ostriker & Peebles 1973, White 
1983b). The initial clouds are in polytropic equilibrium with small internal 
velocities, compared to what they would need to be in dynamical, gravita- 
tional equilibrium. This feature avoids a large initial spin, in accordance 
with the fact that there were no primordial rotational motions in the uni- 
verse (Ozernoy & Chernin 1968; Parijskij 1973; Boynton & Partridge 1973; 
Peebles & Silk 1990). Then, the clouds are sent to approach each other, and 
only after their gravitational interaction takes place, spin will be gained. 

For the simulations, each protogalaxy is constructed using the Plummer 
model given by the potential-density pair (Aarseth et al 1974), 



eter which determines the dimensions of the cloud. Particle velocities are 
chosen everywhere isotropic which gives a system initially in steady state. 
The total energy of the cloud is E — — (37r/64)GM 2 6 _1 . We are using units 
in which G = M = -4E = 1. 

We take two identical 3-D clouds consisting of N = 2 11 particles. The 
initial separation between the clouds is 10, and the velocity of the center of 
masses are V\ = (0.1,0.1,0) and V2 = (—0.2,0.4,0) in our units. The initial 
velocities are given so that the kinetic energy is a fraction of the potential 
energy and we consider a range that is consistent with the observed velocities 
of galaxies in clusters that goes from 50 km/s up to about 1000 km/s. For 
the present investigation we consider that the protogalaxies moves initially 
in the plane (x,y), and the angle of both protogalaxies is the same. More 
general initial conditions will be considered in a future communication. 

A three-dimensional hydrodynamic code based on the TREESPH algo- 
rithm formulated by Hernquist & Katz (1989) was used for the computations 
of this paper. The code combines the method of SPH, developed by Lucy 
(1977) and Gingold & Monaghan (1977), with the hierarchical tree algorithm 
of Hernquist (1987) for the calculation of the gravitational acceleration forces. 





GM 





where G is the gravitational constant, M is the total mass and b is a param- 
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Figure 1: Evolution of the angular momentum without scalar field and with 
scalar field for different values of A 
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The SPH method is a grid-free Lagrangian scheme in which the physical prop- 
erties of the fluid are determined from the properties of a finite number of 
particles. In order to reduce the statistical fluctuations resulting from repre- 
senting a fluid continuum by a finite set of particles, a smoothing procedure 
is employed in which the mean value of a field quantity is estimated from 
its local values by kernel interpolation. Thus, the evolution of particle i is 
determined by solving Euler's equation 

v, (10) 
--Vpi - \v(h m )i + A visC)i , (11) 

where pi and A visc ^ denote, respectively, the gas pressure and the artificial 
viscous acceleration associated with particle i. This quantities are introduced 
because we are considering that the protogalaxies are gaseous. The code 
was modified (Rodrfguez-Meza 2000) to include the effect of the scalar field, 
through the term h 00 given by Eq. (6). 

The simulation of the interaction of two protogalactic models starts with 
the clouds separated by a distance of 10 and on the The selected 

separation is large enough to ensure that tidal effects are important but 
small enough that the calculation is possible in a reasonable computing time. 
Each particle in the initial steady state clouds is given an additional velocity 
(Vi or V2, corresponding to cloud 1 or 2) so that its magnitude is much 
bigger than the internal velocities they have at the equilibrium described 
above. In this way, initial clouds are almost spinless, and the given velocities 
Vi, V 2 imply kinetic energies associated with each cloud. The evolution of 
the intrinsic angular momentum (J/M) with respect to the center of mass 
of each cloud is shown in Fig. 1. Continuous lines indicate cloud 1 and 
dashed lines cloud 2. The first plot is without scalar field, the other plots 
consider values of A of b/16, 6/32, and 6/64, respectively. We observe that 
the intrinsic angular momenta start from their initial values to a constant 
mean value approximately of 0.75 which in the cgs units is of the order of 
10 30 cm 2 /s. This transient stage is slower without scalar field than the ones 
which consider scalar field. The faster transient occurs when A is bigger. In 
Fig. 2 we show plots of phase space v r versus r of the whole system and for 
the same cases as Fig. 1. The scalar field extends the phase space in the v r 
direction. 



dt 
dVi 

dt 
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Figure 2: Phase space v r versus r of the combined system after the transient 
stage without scalar field and with scalar field for different values of A 
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4 Conclusions 



We have made other computations varying the kinetic energy from 4 to 1/16 
times the potential energy. For large values of the kinetic energy the deflec- 
tion is small, but for small values there is a considerable deflection, and in 
some cases we got almost a head-on collision. This is consistent with the 
known fact that the merging probability in an encounter of two clouds is en- 
hanced significantly when the encounter takes place at relatively low speed 
(see for instance Makino & Hut 1997). We found that only close encounters 
and mergings permit the original spinless clouds to gain rotational velocities 
as is observed in typical galaxies nowadays. Similar studies have being done 
(Namboodiri & Kochhar 1991) considering point mass perturbers. In our 
approach the perturber is itself a protogalaxy and therefore the dynamics is 
more complicate, especially in close encounters. We have also found that the 
transient time to spin up the clouds depends on the scalar field. The tran- 
sient stage is faster than the one without scalar field. When the scalar field 
is included faster transients occur for bigger values of A. The phase space 
v r versus r of the combined system is also more extended in the v r direction 
with scalar field than the one without scalar field. 
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